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Abstract 

The performance of the method of angular moments on the Ar s determination from 
analysis of untagged decays B®(t), B° s (t) — > J/^(— > K + K~) is examined 

by using the SIMUB generator. The results of Monte Carlo studies with evaluation 
of measurement errors are presented. The method of angular moments gives stable 
results for the estimate of Ar s and is found to be an efficient and flexible tool for the 
quantitative investigation of the — > J/i/jcf) decay. The statistical error of the ratio 
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1 Introduction 

The study of decays B°(t),B° s (t) -> J/^H / + H K + K~), which is one of the 
gold plated channels for 5-physics studies at the LHC, looks very interesting from the 
physics point of view. It presents several advantages related to the dynamics of these 
decays, characterized by proper-time-dependent angular distributions, which can be 
described in terms of bilinear combinations of transversity amplitudes. Their time 
evolution involves, besides the values of two transversity amplitudes at the proper 
time t — and their relative strong phases, the following fundamental parameters: 
the difference and average value of decay rates of heavy and light mass eigenstates 
of Bg meson, AT S and r s , respectively, their mass difference AM,, and the CP- 
violating weak phase <f>( s \ The angular analysis of the decays B®(t),B~° s (t) — > J/ij){— > 
K + K~) provides complete determination of the transversity amplitudes 
and, in principle, gives the access to all these parameters. 

In the present paper we examine the performance of the angular-moments method 
[1] applied to the angular analysis of untagged decays B®(t), B° s (t) — > J/ip(— > l + l~) <f>(- 
K + K~) for the determination of Ar s . After giving the physics motivation in Section 
2, we describe in the next section the method of angular moments based on weight- 
ing functions introduced in Ref . [1] . For the case of Ar s determination this method 
is properly modified in Section 4. The SIMUB-package [2] for physics simulation of 
5-meson production and decays has been used for Monte Carlo studies. A general in- 
formation about the SIMUB package one can find in Appendix. The dynamics of the 
decays Bg(t),B^ s (t) — > J/ip(^ l + l~)cj)(-^ K + K~) is described by four-dimensional 
probability density functions depending on decay time and three physical angles. The 
algorithms of multidimensional random number generation have been elaborated and 
then implemented in the package SIMUB to provide tools for Monte Carlo simula- 
tion of sequential two-body decays B®(t), B s (t) — > J/i[)(— > l + l~)(f)(-^ K + K~) in 
accordance with theoretical time-dependent angular distributions. This algorithms 
is discussed in Section 5. In Section 6 we present the results of the Monte Carlo stud- 
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ies of the angular analysis of the decays B®(t), B s (t) — > J/ip(—> <f>{— > K + K~) 
with two sets of weighting functions for angular-moments method and concentrate 
on the evaluation of measurement errors and their dependence on statistics. 



2 Phenomenological description of the decays 

B° s (t)X(t) - l + l-)<P{^ K+K-) 

The angular distributions for decays B°(t),B° s (t) -> J/^(-> l + l~) K+R-) are 
governed by spin-angular correlations (see [3]- [6]) and involve three physically deter- 
mined angles. In case of the so-called helicity frame [5] , which is used in the present 
paper, these angles are defined as follows (see Fig. 1): 

• The z-axis is defined to be the direction of 0-particle in the rest frame of the 
B° s . The x-axis is defined as any arbitrary fixed direction in the plane normal 
to the ^-axis. The y-axis is then fixed uniquely via y = z x x (right-handed 
coordinate system). 

• The angles (0;+, Xi+) specify the direction of the / + in the J/ip rest frame while 
(@k+j Xk+) give the direction of K + in the rest frame. Since the orientation 
of the x-axis is a matter of convention, only the difference x — Xi+ ~ Xk+ of 
the two azimuthal angles is physically meaningful. 

In the most general form the angular distribution for the decay B®(t) — > J/ip(^> 
l + l~) > K + K~) in case of a tagged B° s sample can be expressed as 

<VX "' HB " ) - !) £ Oi(t)gi(e l+ , e K+ , x ). (i) 



dcosQi+ dcosQ K + dx dt 32n i , 
Here Oi (i = 1,...,6) are time-dependent bilinear combinations of the transversity 
amplitudes A (t), A\\(t) and A±(t) for the weak transition B°(t) — > J/ip(/> [7] (we 
treat these combinations as observables) : 

O l = \A,{t)\ 2 , 2 = \A ll (t)\\ 3 = \A ± (t)\ 2 , 



4 = Im(^(*HL(f)), 0, = Re(A*(t)A {l (t)), O, = lm(A* (t)A ± (t)) ,(2) 

and the gi are functions of the angles 6;+, x onr y [5]: 

<7i = 2cos 2 Ox+sin 2 9;+ , 

5r 2 = sin 2 6 x +(l - sin 2 6i+cos 2 x) , 

g 3 = sin 2 6x+(l - sin 2 ©i+sin 2 x) , 

gi = — sin 2 6x+sin 2 6;+sin2x , 

#5 = -^sin26;+sin20^+cosx , 
v2 

g 6 = -^=sin2e i+ sin20 K +sinx. (3) 

For the decay B~° s (t) — > J/^(— > Z + Z~) 0(— > K + K~) in case of a tagged B° s sample the 
angular distribution is given by 

(/ ' v "" uo - 9 T,Vi(t) gi (e l+ ,e K+ ,x) (4) 



dcosQi+dcosQ K +dxdt 32ir ~{ 
with the same angular functions g t and 

1 = |A(t)| 2 , 2 = \A ll (t)\\ 3 = \A ± (t)\\ 

4 = lm(A^t)A ± (t)), 5 = Re(A*(t)A l{ (t)), % = lm(A* (t)A ± (t)) ,(5) 
where Ao(t), and are the transversity amplitudes for the transition 

The time dependence of the transversity amplitudes for the transitions ,B s (t) 
J/ift (f) is not of purely exponential form due to the presence of i?!? — 5° mixing. This 
mixing arises due to either a mass difference or a decay-width difference between the 
mass eigenstates of the (B^ — B° s ) system. The time evolution of the state \B®(t)) 
of an initially, i.e. at time t — 0, present B® meson can be described in general form 
as follows: 

\B° s (t)) = g + (t)\B° s ) + g-(t)\B° s } , g + (t = 0) = 1 , g_(t = 0) = 0, 
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i.e., the state \B®(t)) at time t is a mixture of the flavor states \B®) and \~B^ S ) with 
probabilities defined by the functions g+(t) and g~(t). In analogous way, the time 
evolution of the state \~B° s (t)) of an initially present B~° s meson is described by the 
relation 

= 9 + (t)W) + 9-(t)0l) , 9 + (t = 0)=0, M* = 0) = 1. 
Diagonalization of the full Hamiltonian (see [8] for more details) gives 

9+(t) = \(f^ Lt + e-* Bt ) , g-(t) = |(e-^* - e"*-') , 

S+(t) =g_(t)/a 2 , g_(t) = g + (t) . (6) 

Here Hl/h = ML/H — (i/i)TL/H are eigenvalues of the full Hamiltonian corresponding 
to the masses and total widths of "light" and "heavy" eigenstates I-Bl/jj), and a is 
a phase factor defining the CP transformation of flavor eigenstates of the neutral 
£? s -meson system: CP\B° S ) = a\B° s ). In the case \a\ 7^ 1 the probability for B® to 
oscillate to a B° s is not equal to the probability of a B° s to oscillate to a Such an 
asymmetry in mixing is often referred to as indirect CP violation, which is negligibly 
small in case of the neutral B-meson system. 

The time evolution of the transversity amplitudes Af(t) (/ = 0, ||, _L) is given by 
the equations 
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A f (t) = A f (0) g + (t) + g_tt)-^tf , A f {t) = A f (0) g + (t) + // H)—C f 



Vcp 01 
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V C p a 



(7) 

Here r) CP are eigenvalues of CP-operator acting on the transversity components of 
the final state which are eigenstates of CP-operator: 

CP\J/^<t>) f = v f C p\J/^<t>)f , (/ = 0, ||, J.) , 
Vcp = 1 > Vcp = 1 > Vcp = ~ 1 > 
and £^ is the CP- violating weak phase [9]: 



^ = e-^ S) , 0« = 2[arg(y t ;y t6 ) - arg(^)] = -25 7 , 



where 5 is the complex phase in the standard parameterization of the CKM matrix 
elements V^- (i G {u, c, t}, j G {d,s,b}), and 7 is the third angle of the unitarity 
triangle. 

The phase <f>^ is very small and vanishes at leading order in the Wolfenstein 
expansion. Taking into account higher-order terms in the Wolfenstein parameter 
A = sin^c = 0.22 gives a non- vanishing result [10]: 

0(s) = _2X 2 V = -2A 2 i4sin 7 . 

Here 



Rb = 



1 


\v ub \ 


A 


\Vcb\ 



is constrained by present experimental data as R b = 0.36 ± 0.08 [11]. Using the 
estimate 7 = (59 ± 13)° [12], the following constrain can be obtained for the phase 
#: 

0« = -0.03 ± 0.01 . (8) 

According to Eq. (7) at time t — 0, the transversity amplitudes of B® , B° s — > 
J/ip(f> decays depend on the same observables |A (0)|, |A||(0)|, |Aj_(0)| and on the 
two CP-conserving strong phases, S 1 = ar g[A^(0)A±(0)} and S 2 = arg[AQ(0)A±(0)}. 
Time-reversal invariance of strong interactions forces the form factors parameterizing 
quark currents to be all relatively real and, consequently, naive factorization leads 
to the following common properties of the observables: 

ImL4*(0)A ± (0)] = 0, Im[^(0)A ± (0)] = 0, Re[A* (0)A {{ (0)} = ±\A (0)A {{ (0)\ . 

Moreover, in the absence of strong final-state interactions, 5i = n and 5 2 = 0. 

In the framework of the effective Hamiltonian approach the two body decays, 
both — > J/ip(f) and B® — > J/ipK*, correspond to the transitions b — > sec with 
topologies of color-suppressed spectator diagrams shown in Fig. 2. Factorizing the 
hadronic matrix elements of the four-quark operators of the effective Hamiltonian 
into hadronic matrix elements of quark currents, the transversity amplitudes |A (0)|, 
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|A||(0)|, \A ± (0)\ of decays B^,W q -> J/ipV((q,V) G {(s, 0), (d, K*)}) can be ex- 
pressed in terms of effective Wilson coefficient functions, constants of J ftp decay, 
and form factors of transitions B q — > V induced by quark currents [1]. In Table 1 
we collect the predictions of Ref. [1] for the transversity amplitudes of B® — > J/ip(f> 
(B% — » J/ipK*) calculated with B — > X* form factors given by different models 
[13, 14, 15]. The B ^ K* form factors can be related to the B — > case by us- 
ing SU(3) flavor symmetry. The most precise polarization measurements performed 
recently in decays S — > J/ipK*: 

|A (0)| 2 = 0.60 ±0.04, |A ± (0)| 2 = 0.16 ±0.03 (BaBar [16]) , 
|A (0)| 2 = 0.62 ±0.04, |A ± (0)| 2 = 0.19 ±0.04 (Belle [17]), 

confirm the predictions based on the model [15]. 



3 Angular- moments method 

The angular distributions for decays B®(t), B° s (t) — > J/ip{— > Z + Z~) 0(— > K + K~) in 
case of tagged 5° and .B s (i) samples (see Eqs. (1) and (4), respectively) as well as 
in case of the untagged sample can be expressed in the most general form in terms 
of observables bi(t): 

f (e l+ ,e K+ ,x;t) = ^- E 6 <(*)ft(©i+>©*r+>x)- (9) 

32tt i=1 

The explicit time dependence of observables is given by the following relations: 



h(t) = 


|4>(0)| 2 G L (f), 


hit) = 


|A„(0)| 2 G L (f), 


hit) = 


\A ± iO)\*G H it), 


hit) = 


1^1(0)11^(0)1^), 


hit) = 


\Ao(p)\\M Q )\ G L(t) COs(5 2 -5l), 


hit) = 


\A iO)\\A ± iO)\Z 2 it), 
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where we have used the general compact notations: 

Gl/hH) = l - [(1 ± cos0W)e- r ^ + (1 =F cos^e"™ 
Zi, 2 (t) = l(e- THt ~ e- r ^)cos5 li2 sin0^ 



- for observables hi = (Oi + Oj)/2 in case of the untagged sample with equal initial 
numbers of and T3° s , while 

GffV&\t) = G L/H (t) ± e- r ^sin(AMt) sm^ , 



r L/H 

z{f )nB "\t) = Z 1>2 {t) ± e~ rs * [sin5 li2 cos(AMt) - cos5i, 2 sin(AMt)sin^ 

- for observables 6- Bs ' ) = 0, and 6- Bs ' ) = Oj in case of tagged B° s and B° s (t) samples, 
respectively, with T s = (r L + Y H )/2. It is easy to see that both in the tagged and 
untagged case we have 

G L/ H(t)\^ =0 = e- r ^. 

According to Ref. [1], the observables bi(t) can be extracted from distribution 
function (9) by means of weighting functions Wi(Qi+, Qr+, x) f° r eacn i sucn that 
9 f 

— J dcos<d l+ dcosQ K +dx Wi(<d l+ ,<d K+ ,x) 9j(&i+, x) = $ij , ( n ) 
projecting out the desired observable alone: 

hit) = J dcosQi+dcosQ K +dx /o(Oj+, @k+, X 5 ^(9,+ , 9^+, x) • (12) 
The angular-distribution function (9) obeys the condition 

L(t) = J dcosQ l+ dcosQ K+ d X / (e,+ , 9^+, x ; *) = &i(*) + fc 2 (0 + & 3 (*) • (13) 

For decays B — > J/ip(— > l + l~)<fi(-^ K + K~), the explicit expressions of weight- 
ing functions, given in Table 5 of Ref. [1] for physically meaningful angles in the 
transversity frame, get the following form (Set A) after transformation into the he- 
licity frame: 

W { A ^ = 2-5 cos 2 9 /+ , 



(A) 2 2 

w 2 = 2 — 5 sin 6/+ cos x > 

w ^ = 2-5 sin 2 9 /+ sin 2 x , 

(A) 5 . 2 
u;} = -- sin e x +sm2x, 

(A) 25 

= sin 2<d K + sin 29/+ cos x , 

25 

Wg 4 '* = — j= sin 26^+ sin 26/+ sin x ■ (14) 

The expressions of Eq. (14) are not unique and there are many legitimate choices 
of weighting functions. A particular set can be derived by linear combination of 
angular functions g t (see [1] for more discussions): 

6 

w i (Q l +,Q K +, x) = Yl \j9j{®i+,®K+, x) , (15) 

where the 36 unknown coefficients \j are solutions of 36 equations 
9 



\j J dcosQi+ dcosQ K + dx g j (Qi+,Q K +, x) 9k(®i+,®K+, X) = $ik • (16) 



32vr j=1 

The weighting functions (set B) corresponding to the linear combination of the an- 
gular functions (3) are given by 

W W = -L[28cos 2 9 x +sin 2 9/+ - 3 sin 2 9^+(l + cos 2 9/+)] , 

w ( 2 B) = -i[4cos 2 e x +sin 2 0;+ - 29 sin 2 9^+ (1 - sin 2 9,+ cos 2 x) 
8 

+ 21 sin 2 Q K +(1- sin 2 Q l+ sin 2 x)\ , 

w ( 3 B) = -i[4cos 2 9 x +sin 2 9 i+ +21sin 2 9 x +(l -sin 2 9 ; +cos 2 x) 
8 

- 29 sin 2 e K +(l - sin 2 Q l+ sin 2 x)} , 
25 

= sin 2 9^+ sin 2 9 ; + sin 2% , 

8 



wi ' — w 



5 ) 

4 B) = wi A K (17) 

For a limited number of experimental events N in the time bin around the fixed 
value of the proper time t, distributed according to the angular function (9), it is 



convenient to introduce the normalized observables 

1,(1)= b^t) /L{t) (18) 

with normalization factor L(t) given by Eq. (13). Then, as it follows from the 
Eq. (12), the observables bi{t) (18) are measured experimentally by 

1 N 

JV 3=1 

with summation over events in a time bin around t. Here w\ = Wi(Q 3 l+ , ® J K +, X*)-> 
where 6^+, 0-^+ and x J are angles measured in the j-th event. The statistical 
measurement error of the observable (19) can be estimated as 



N' 



N 

, z(r p) -w) 2 , 

V=i 



with summation over all events in the same time bin. 



4 Time- integrated observables 

For data analysis it is rather convenient to use the time-integrated observables [18] 
defined as 

1 f T o r 

k(T ) = / dt / dcos<di+dcos<d K +dxWi(<di+,<d K +,x) /o(©/+, X\ *) 

L(T) jo j 

(20) 

with argument T < T, where T is the maximal value of the 5-meson proper time 
measured for the sample of events being used, and L(T) is a new normalization 
factor, which has the form: 

L(T) = J L(t) — J dt J dcosQi+ dcosQ K + d\ /o(cosOj+, cos6^:+, x\ t) — 

= (|A (0)| 2 + |An(0)| 2 )G L (T) + |A ± (0)| 2 ^(T), (21) 

where, in the compact notations used in Eq. (10), 

T 



G L/H (T) = jf dtG L/H (t). 



The following normalization condition is valid for the observables (20): bi+b 2 +b 3 = 1. 
For a limited number of experimental events N(T), measured in the proper time 
region t G [0,T], Eq. (20) reduces to 



r(exp) 



(To) 



N(T ) 

E 



(22) 



N(T) . , 

with summation over all events N(T ) in the time interval t G [0, T ] for T < T. In 
case of the untagged sample we have 



Gl/h(T) 



,-rvT 



+ (1 =F cos0^) 



(23) 



and 



Z(T) = 



cos^i i2 sin0 



^^i, 2 (T) 



e-r ff T _ 1)AV _ ^ _ 1)/r 



Y,T 



For the untagged sample the explicit form of time-integrated normalized observ- 
ables (20) in terms of the functions Gl/h{T) and Z(T) is given by 

b 1 (T ) = \A (0)\ 2 G L (T )/L(T), 

b 2 (T ) = \A ll (0)\ 2 G L (T )/L(T), 

h(T ) = \A ± (0)\ 2 G H (T )/L(T), 

b 4 (T ) = |A,|(0)| \A ± (0)\ Z(T o )cos5 1 sin0( s )/L(T) , 

6 5 (T ) = \A (0)\ |A,|(0)| G L (T ) cos(5 2 - <Ji)/L(T) , 

6 6 (T ) = |4o(0)| |A ± (0)| Z(T ) cos5 2 sin0( s V^(^) • (24) 

In the Standard Model (SM) sin0^ xs and the observables b4 : e(T ) are vanishing. 
In case of a new physics signal the values of sin0^ and 64,6 (T ) can be sizable, 
however. 

The following relations are valid for the observables (24): 



h(T ) = cos5 1 sin0«Z(To 



b 2 (T )~b 3 (T ) 
\G L (T )G H (T ) 
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b 5 (T ) = cos(5 2 - S,) V&i(T ) b 2 (T ) , 
b 6 (T ) = cos<5 2 sin0^Z(T o ) 
If we introduce the function 



61 (Tp) 6 3 (Tp) 
^G l (T )Gh(T ) 



7(T) = G H (T)/G L (T) , (25) 

then, the values of initial transversity amplitudes at t — and the strong-phase 
difference (5 2 — 5i) are determined from the observables bi(T) = 6j(T = T ) by 



|4>(0)|< 



6!(T) + 6 2 (T) + 6 3 (T)/7(T) 



l^||(0)| 2 =~ 



6i(T) + 62(T)+6 3 (T)/7(T)' 

l^(0)| 2 



6 3 (T)/7(T) 



6i(T)+6 2 (T)+6 3 (T)/7(T) 

cos(5 2 - 5 X ) = &s(T) , (26) 
A CO 62 CO 

where we consider the initial amplitudes normalized as |A (0)| 2 +|A||(0)| 2 +|A_|_(0)| 2 = 
1. We have also: 

^ cmSl , = . M3 . (27) 

Ai(T)5a(T) ^ 

For extraction of the _E>°-width difference Ar s = r# — T L from experimental data 

it is convenient to use a special set of the time-integrated normalized observables 

[18]: 

6»(T ) = / dt / dcos0i+dcos0 K +dx^(@i+ ) @K+,x) e r '* / o (0j+, ©x+, x; *) , 

L(T) Jo j 

(28) 

where T' is some arbitrary initial approximation of the 5°-meson total decay width. 
These observables can be extracted from the experimental events N(T), measured 
in the proper time region t G [0,T], by using the formula 



1 



N(T ) 



r P) (To) = ^y £ WU (29) 
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where W- = e r ' t3 wj, and summation is performed over all events N(T ) in the time 
interval P E [0,T ]. 

For the untagged sample, the explicit expressions for the time-integrated observ- 
ables (28) can be easily obtained by replacing bi, Gl/h and Z in the expressions of 
Eq. (24) by hi, Gl/h and Z, respectively, (with the same normalization factor (21)) 
after introducing the following notations 



G L/H {T) = J* dte r,t G L/H (t) 



AT L T/2 _ i ~AV H T/2 _ i 

= (lioos^W)— ^ (l=Foos^W) ^- , (30) 

ZiT) ee - - f dt^Z^T) 
cos0i i2 sm<pc Jo 

1 _ e Ar L T/2 x _ e -AY H T/2 

AT\ + ATV ' 

where AT l /h are auxiliary parameters given by 

Ar L = 2(r - v L ) , Ar H = -2(r' - v H ) . (31) 

Eq. (26) is also valid after such a replacement. 



5 Algorithms of multidimensional random num- 
ber generation 

The time-dependent angular distribution function (9) is used in the SIMUB gener- 
ator [2] (see also Appendix) as density of the probability function for Monte Carlo 
simulation of the vertex and kinematics of the final-state particles in case of de- 
cays B° s (t),B° s (t) -> J/V>H l + l~)4>{-^ K + K~). The variables of the function 
/o(©z+, 6jf+, Xi t) can n °t be factorized and randomly generated in an indepen- 
dent way. Nevertheless, the random generation of cos0/+, cosO/^+, x an d t can 
be performed either simultaneously according to the distribution function (9) by 
using four-dimensional random generator or successively, one after another, by us- 
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ing single-dimensional random number generators with accordance to distribution 
functions obtained by successive integration of the function (9) over its variables. 

Let as consider the latter approach in the case of sequential random generation 
of the variables t, X , cos 6/+ and cos0^+. The following three distribution functions 
are used in this case: 

r+l 

/i(©j+,X;*) = J dcosQ K + fo(Oi+,Q K +,x;t) , 

f 2 ( X ; t) = £ dcosB l+ h(e l+ , X \ t) , (32) 

and 

f 3 (t)= / d X f 2 ( X ;t) =b 1 (t)+b 2 (t) + b 3 (t). (33) 
Jo 

The functions of Eqs. (32) and (33) present the most important experimentally ob- 
servable distributions. In case of sequential decays B° s (t),jT s (t) -> J/^(-> 

K + K~) the explicit form of the distribution functions (32) are given by 

4 4 
/i(cos0i+, x\ t) = ^h(t) sin 2 9 i+ + -62(f) (1 - sin 2 /+ cos 2 x) 

+ ^b 3 (t) (1 - sin 2 /+ sin 2 x) - C 4 (t) sin 2 /+ sin2x , 

f 2 ( X ,t) = b l (t) + h 2 (t)(3-2cos 2 X ) 

+ h 3 (t) (3 - 2sin 2 X ) - b A (t) sin2 X . (34) 

The procedure of random generation of the variables t, X , cos 0^+ and cos 0^+ is 
as follows: 

• First, the proper time t is randomly generated according to the distribution 
function /3(f). 

• Second, the angle X is randomly generated according to the single-dimensional 
distribution given by function f 2 ( X ;t) with t being fixed to be equal to the 
time, generated at the first step. 

• Then, the value cos Z + is generated according to distribution function f 1 (0 ; + , X ; t) 
with values of t and X fixed to be equal to their values generated at the previous 
two steps. 
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• Finally, the value of cos &k+ is generated according to the single-dimensional 
distribution given by function /o(@z+, ®k+, x; £) with a properly fixed values 
of t, x an d cos 9;+. 

Two Monte Carlo methods of random numbers generation are implemented in 
the C++ class T_VertexB_VllVpp of the BB_dec program of the SIMUB package [2]. 

The first method is based on filling of the single-dimensional array f 4_Integ [n4_cells] 
of real numbers which represents numerically the four-dimensional distribution func- 
tion /o(©z+, X; t). The array f4_Integ contains n4_cells = iV(cos0;+) x 
iV(cos©^+) x N(x) x N(t) elements, where iV(cos0 i+ ), iV(cos0 x +), N(x) and N(t) 
are the numbers of bins (generator resolutions) for corresponding variables randomly 
generated in a four-dimensional volume V G {cos©;+, cosO^+, x, t}. A fast gener- 
ation of these variables is performed by the function GetRandom4 by using the al- 
gorithm which is analogous to that was realized in the ROOT class TF3 [19] for 
generation of three random numbers distributed according to a three-dimensional 
probability function. In this case the array f 4_Integ is filled in constructor of the 
class T_VertexB_VllVpp where a large memory space for the array is reserved during 
all time of existence of this class sample. 

The second method involves sequential generation of random numbers according 
to approach based on usage of the single-dimensional distribution functions given 
by Eqs. (33) and (34). The corresponding algorithm of generation of the vari- 
ables t, x? cosO/+ and cosQ K + is described above. In this case the constructor 
T_VertexB_VllVpp does not fill any arrays and, therefore, it does not reserve a large 
memory space. 

Both methods demonstrated the identical results with reasonable CPU time and 
memory used. 

The .B-decay dynamics described by angular correlations (9) should be included 
without fail into Monte Carlo generators developed for physics studies. The main mo- 
tivation for developing the new dedicated package SIMUB was that already existing 
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generators do not take into account all theoretical refinements which are of a great im- 
portance for Monte Carlo studies of 5-decay dynamics. In particular, in the genera- 
tors PYTHIA [20] and QQ [21] the time- dependent spin-angular correlations between 
the final-state particles in the decays Bg(t),B^ s (t) — > J/ij}{— > l + l~) (f>{— > K + K~) are 
not reproduced in a proper way 1 . 

The time and angular distributions in the helicity frame for decay B®(t), B s (t) — > 
J/^(-> A* + A*~) K + K~), generated by SIMUB, PYTHIA and QQ packages, are 
compared in Fig. 3. In case of PYTHIA usage, Fig. 3 shows the uniform distribu- 
tions for angular variables cos 0x+, cos 6^+ and x because of lack of time-dependent 
angular correlations. Due to this simple reason PYTHIA can not be used for Monte 
Carlo studies of dynamics of sequential two body decays of B mesons in the channels 
of the type B — > Vi(— > At + At~) V 2 (^ P + P ) with intermediate vector mesons V\ and 
V 2 . Unfortunately, another well known package QQ turns out to be also not suitable 
for study of decays of this type because of lack of azimuthal-angle x correlations. 

6 Monte Carlo studies 

For Monte Carlo studies of the estimation of physical parameters by applying the 
angular-moments method, untagged samples of events of the decays B®(t), B s (t) — > 
J/ip(fi have been generated by using the package SIMUB with various sets of the 
input values of initial amplitudes |A)(0)| and |Aj_(0)| and AT S . Other parameters 
are fixed as follows: 

S x = 7T , 5 2 = , T s = 1/t, = 2.278 [mm/c]" 1 , <p{ s) = 0.04 . 

The value of r s used corresponds to the lifetime r s = 1.464 ps [12] while the CP- 
violating weak phase <f>^ was fixed as the upper limit of the constrain (8). The value 

of Ar s is expected to be negative in the SM. The combined experimental result 
1 In the latest version of the generator EvtGen [22] this mode can be also simulated including 
the full time-dependent spin-angular correlations. 
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for | Ar s | /T s is not precise: | AT S | /T s < 0.52 at 95% CL [12]. In the approxima- 
tion of the equal B° s and B® lifetimes, the |Ar s | extraction can be improved [12]: 
| Ar s |/r s < 0.31 at 95% CL. A set of the untagged-event samples has been generated 
with AT S /T S G [—0.3, —0.01] to study the influence of AT S value on the estimation 
of Bg(t) — > J/ip(f> decay parameters from data analysis. 

The values of the time integrated observables Vf xp \T G ), defined by Eq. (20), can 
be extracted from data according to Eq. (22) by summation of weighting functions 
for each event. The statistical error of 6j(T ) is defined by 



(Sh) 



\(stat) 



N(T)\ 



N(T ) 
3=1 



(35) 



while a systematic error due to limited precision of angular measurements can be 
estimated as 

{5bi) {sys) 



\ N(T) 



(36) 



Here 



dwi 



1 2 



(9cos©i+ 



AfcosG 



+ 



-A(cos9 x +) 



dx 



A(x) 



1 2 



.0 COS <S>K+ 

In a similar way, the values of the observables b[ exp \T ), defined by Eq. (28), 
can be extracted from the data according to Eq. (29). The formulae for statistical 
and systematic errors for bt P) (T ) can be obtained by replacement of w\ to W/ in 
Eq. (35) and the following redefinition of A\ in Eq. (36): 



d cos 6z+ 



A(cos 



-A(cos Qr+) 



+ 



A(x) 



+ 



dt 



A(t) 



_<9cos Qk+ 

Eq. (36) can be applied to estimate the systematic errors related both to the mea- 
surement precision of the detector and to the limited resolution of the Monte Carlo 
generator. In the SIMUB generator, for each variable V G {cos©;+, cos©^+, t} 
randomly generated for decays B®(t) — > J/^(— > Z + /~) 0(— > K + K~), the number of 
bins in the region [V min , V max ] was set as N gen = 50 000. The generation precision for 
the variable V is defined as A(V) = (V max — V min )/N 9en and systematic errors (36) 
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are proportional to (jV 96 ™)^ 1 / 2 . The 5-meson proper time was generated within the 
interval t £ [0, T — 2mm/c] which includes 99.3% of all £?-decays. We have used 
samples with a maximum of 100 000 events of the decay B° s — > J/V0 because a 
statistics of about 80 000 events is expected to be obtained per year at the CMS 
detector at the LHC low luminosity under realistic triggering conditions [23]. 

Table 2 shows the values of the observables bf xp) {T) = bi exp) (T = T) extracted 
from the Monte Carlo data by applying the sets A and B of weighting functions, 
given by Eqs. (14) and (17), respectively. Various theoretical models for estimation 
of the transversity amplitudes |A)(0)| and |Aj_(0)| (see Table 1) have been considered 
to fix these parameters in the SIMUB generator. It can be seen from Table 2 that the 
choice of N gen = 50 000 provides negligibly small systematic errors for the observables 
as compared with the statistical ones. Moreover, both errors slightly depend on the 
values of the observables. For observables obtained by using the set-B weighting 
functions, the statistical errors are significantly smaller than in case of the set-A 
weighting functions. We should also note that even with the statistics of 100 000 
events, the values of observables (T) and - as consequence of Eq. (27) - the 
combination cos5i )2 sin0 cannot be extracted from the data if the CP-violating weak 
phase (f>W is small according to the SM expectation (8). In this case, these parameters 
can be estimated only by using a statistics which is not less than 3 x 10 9 B®(t) — > 
J/ip(/> decays. 

Analysis of the same Monte Carlo data leads to similar conclusions concern- 
ing the behavior of statistical and systematic errors for the observables bf xp \T) = 
b\ exp \T = T). To illustrate the performance of our method in this case, only the 
results obtained for transversity amplitudes, corresponding to the Cheng's model 
[15], are shown in Table 3. 

Fig. 4 shows the dependence of the observables h(T) and 

bm= l ~ T e ~^ T b t (T) (i = 1,2,3) 
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on the value of the ratio AT S /T S . For Ar s = we have 

6i A s(T)|at.=o = ^asCOIam = |4>,||,±(0)| 2 . 

The observables bi(T) slightly depend on Ar s . The rather strong dependence of the 
observables h(T) on the decay width difference Ar s , shown in Fig. 4, can be used 
for extraction of this parameter from the data analysis as it will be discussed below. 
Under the assumption <f>^ = 0, we have from Eq. (30): 

~ m e ±AF L/H T/2 _ j 

G?L(T) = ±2- — . 

L/HK ' at lih 

Therefore, the values of the auxiliary parameters AT L / H , defined by Eq. (31), can 
be determined separately by using the ratios of observables b\ exp \T) /b\ exp \T Q ) , ex- 
tracted from the data analysis, and solving numerically the equations which arise 
from one of the following relations: 

km/kiTo) = G L (T)/G L (T ) (i = 1, 2, 5) (37) 

- to determine Ar^, and the relation 

b 3 (T)/b 3 (T ) = G H (T)/G H (T ) (38) 

- to determine ATh- Then, the decay-width parameters T s , AT S and Tl/h can be 
determined via V and Ar L /^ as 

Ar L - AT H AT L + AT H r^ AT ^ 
1 s = 1 ^ , Al s = , 1 L/H = 1 =p ■ (39) 

So, by using some reasonable approximation for V as a starting point for the data 
analysis, the experimental value of Y s can be essentially improved simultaneously 
with determination of Ar s . The statistical error of T s determination is expected to 
be twice smaller than for Ar s determination. 

The direct numerical calculations have shown that the difference between the 
values of observables &«(T) (i = 1,2,3,5), calculated with 0^ — o and = 0.04, 
does not exceed 0.01%. Even in case of statistics of 100 000 events this difference is 
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negligibly small as compared with statistical errors for these observables (see Table 3). 
Therefore, the assumption = is a good approximation for T S) AT S and Y L / H 
determination by the method described above. 

Table 4 shows the results of determination of the decay-width parameters after 
applying the described procedure to the Monte Carlo data. The sample of 100 000 
events generated in case of Cheng's model with Ar s /r s = —0.15 has been used. Both 
sets A and B of weighting functions have been applied to extract the observables 
b\ exp \ The value of T', which is treated as some arbitrary initial approximation for 
the total decay width of _B°-meson, was fixed as V = 1.05 T s , i.e. it was shifted by 
5% relative to the "true" value of T s fixed in the Monte Carlo generator SIMUB. 
The value of To = 0.1 T was chosen as it provides the minimal statistical errors 
to determine the ratios b[ exp \T) /^^^(Tq) . In Table 4 we present the result for 
Ar L obtained from the ratio b[ exp \T)/b[ exp \T ) only, which gives the best precision. 
Table 4 shows that the set-B weighting functions give more precise and stable results 
than the set-A functions. 

To improve the precision of Ar s determination, the same procedure should be 
repeated with T 1 fixed to be equal to the value of T s determined at the first step. Be- 
cause of AT S = ATl = ATh in case of T' = T s , the value of AT S is defined at the sec- 
ond step to be equal to the value of ATl determined from the ratio b^ xp \T) /Ui XP \Tq) 
using Eq. (37). Using the values of AT S from Table 4 as an input value of V at the 
second step, we have obtained finally the following results (to be compared with the 
input value Ar s = —0.3418 set in the SIMUB generator): 

AT ex P = _ .33o ± 0.057 ± 0.004 (set A) , 
AT ex P = _ .338 ± 0.034 ± 0.002 (set B) . 

This way one can reduce not only the statistical error but also essentially improve 
the stability of the Ar s result even in case of using the set-A weighting functions. 

Table 5 shows the statistical errors of AT S /T S determination by the described 
approach applied to different statistics of Monte Carlo events generated with various 
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"true" values of Ar s . The lack of numbers in the table corresponds to cases when 
the approach is not able to give a certain result for Ar s . The use of set-B weighting 
functions gives more stable results even in case of too small statistics and values 
of AT S , for which the same approach does not work with set-A functions. The 
sensitivity of the method is measured by the statistical error of Ar s /r s , which only 
slightly depends on the value of this ratio and is proportional to 1/y/N, where N 
is the number of events. In particular, for a statistics 100 000 events, the statistical 
error is about 0.015, while for 1000 events - about 0.15. 

The value of Ar s can be determined similarly by using the ratios Vf xp \T) /b\ exp \T ) 
or bf xp \T)/Vf xp \T), extracted from the data analysis with V = T s , and solving the 
equations arising from the relations 

b t (T)/b t (T ) = G L (T)/G L (T ) (i = 1, 2, 5) , b 3 (T)/b 3 (T ) = G H (T)/G H (T ) 

or 

k(T)/k(T) = G L (T)/G L (T) (i = 1, 2, 5) , b 3 (T)/b 3 (T) = G H (T)/G H (T) . 

But in both these cases the precision of AT S determination turns out to be worse 
than in the approach based on the ratios bf xp \T) /bi exp \T ) because of the weak 
Ar s -dependence of the bi(T) observables. 

The initial transversity amplitudes and strong-phase difference can be recalcu- 
lated from the values of observables Vf xp \T) according to Eq. (26). The results of 
such determination of the parameters |Aj(0)| 2 (/ = 0, 1 1, _L) and cos(5 2 — ^i) are shown 
in Table 6 for different statistics. We have used the Monte Carlo sample generated 
with the theoretical values of the amplitudes |A (0)| and |Aj_(0)| corresponding to 
Cheng's model [15]. To extract the observables b\ exp \T), the set B of the weighting 
function has been applied to Monte Carlo data. To estimate the statistical errors for 
parameters |Aj(0)| 2 (/ = 0, ||,-L) and cos(5 2 — <5i), the standard error-propagation 
method has been applied to the statistical errors of the observables b[ exp \T), taking 
into account the correlation between pairs of different observables. The systematic 
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errors of the observables related to the limited generator resolution are neglected. 
The total errors for parameters |Aj(0)| 2 (/ = 0, ||,_L) should also include the addi- 
tional uncertainty related to the error of calculation of 7(T) caused by the error of 
Ar s (see definition of j(T) in Eq. (25) and Eq. (23)). In Table 6 we also show these 
errors calculated by assuming Ar s = — 0.15r s (see Table 5) 



7 Conclusion 

For the decay B° s — > J/ip<f> in the framework of the method of angular moments a 
non-fit scheme for separate estimation of the parameters Ar s , T s and |A/(0)| 2 (/ = 
0, ||, -L) has been proposed, based on analysis of an untagged sample, and studied 
by Monte Carlo method. A strong dependence of statistical measurement errors on 
the choice of the weighting functions has been demonstrated. The statistical error of 
the ratio Ar s /r s for values in the interval [0.03, 0.3] was found to be independent of 
the central value and amounts to about 0.015 for 10 5 events. The method of angular 
moments gives stable results for the estimate of Ar s and is found to be an efficient 
and flexible tool for the quantitative investigation of the B° s — > J/ip(/> decay. 
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A General information on the SIMUB package 

The generator SIMUB was developed for Monte Carlo simulation of 5-meson pro- 
duction and decays at the CMS detector. The SIMUB package provides two regimes 
of .B-decay generation: 

• simulation of decay after fragmentation according to angular distribution and 
time dependence (in case of a neutral 5-meson) governed by spin-angular cor- 
relations, time evolution and other theoretical refinements - dynamical mode, 

• simulation of decay by PYTHIA into phase space - PYTHIA kinematical mode. 
The package is kept under the directory SIMUB which has the following structure: 

• bb_gen - routines needed to generate events with 66-pairs at parton level (FOR- 
TRAN, PYTHIA, HBOOK); 

• bb_frg - routines performing string fragmentation, generation of 5-mesons, 
and decays of particles in the PYTHIA kinematical mode (FORTRAN, PYTHIA, 
HBOOK); storing the results into standard HEPEVT or PYJETS Ntuples for 
further usage in the CMS simulation or for analysis; 

• BB_dec - routines performing 5-decays in dynamical mode (C++, ROOT); 

• include - a collection of common blocks for the routines bb_gen, bb_f rg and 
BB_dec; 

• lib - the PYTHIA source codes and object files; 

• doc - documentation. 

The data flow between the three main parts (bb_gen, bb_f rg, and BB_dec) of the 
package SIMUB is shown in Fig. 5 and organized as follows: 

• At the first stage the events containing 66-pairs are generated by the program 
bb_gen, stopped before fragmentation and written to an intermediate Ntuple 
to store events at parton level. 
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• At the next step the program bb_frg reads bb-events from the intermediate 
Ntuple and performs the string fragmentation and the decays of all particles in 
the PYTHIA kinematical mode with exception of _B-mesons selected according 
to user directives to be decayed in dynamical mode. 

• Then, the information about the selected 5-mesons is used by the program 
BB_dec to perform their decays according to dynamical modes. 

• Finally, the information about the selected S-mesons and products of their 
decays is transfered to the subprogram bb_f rg to be added to the information 
about other particles and stored in a final output Ntuple. 

The package SIMUB is installed and tested on Linux (RedHat 6.x, 7.x) platforms. 
The package and documentation for user is found on the SIMUB Package Home Page: 

http://cmsdoc.cern.ch/~shulga/SIMUB/SIMUB.html. 

To install the program, the file SIMUB. tar. gz should be copied from the SIMUB 
Package Home Page to the user directory and unpacked to obtain the directory 
SIMUB. 

The current version of the program SIMUB is adopted to the usage of the gen- 
erator PYTHIA [20] (version 6.215), CERNLIB-2002 [24], and ROOT package [19] 
(version 3.05/07). The standard CERN Program Library should be installed at 
a Linux machine according to the record in the files bb_gen/mak/Makef ile and 
bb_f rg/mak/Makef ile: 

CERNLIB := 'cernlib pawlib graflib mathlib packlib kernlib' 
or the record in the file BB_dec/mak/Makef ile: 
SYSLIBS = -L/cer/new/lib/ -lpacklib . . . 

The user has to generate the object file of the PYTHIA generator in the directory 
lib and to install the ROOT package. Detailed instructions on the installation, 
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compilation and running of the package SIMUB are given in the documents "User 
and Developer Guide" and "Quick Start". The postscript files of these documents 
are placed in the directory SIMUB/doc and also found on the SIMUB Package Home 
Page. 
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Table 1: Predictions for — > J/ipcj) (in brackets - for B% — > J/ipK*) observables 
obtained in Ref. [1] for various model estimates of the B — > X* form factors [13, 14, 
15] (the normalization condition |A (0)| 2 + |A||(0)| 2 + |A^(0)| 2 = 1 is implied) 



Observable 


BSW [13] 


Soares [14] 


Cheng [15] 




^o(0)| 2 


0.55 (0.57) 


0.41 (0.42) 


0.54 (0.56) 




^(0)| 2 


0.09 (0.09) 


0.32 (0.33) 


0.16 (0.16) 
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Table 2: Comparison of the observables b\ exp '(T), extracted from Monte Carlo data, 
with their values bf h \T) corresponding to various theoretical models for |A)(0)| and 
|A_|_(0)|. A sample of 100 000 decay events generated with Ar s /r s = —0.15 was used. 
The first errors are statistical (see Eq. (35)) while the second errors correspond to 
the systematic uncertainties only from limited angular precision (see Eq. (36)) 



a) BSW model [13]: 



i 


bf h) (T) 


b\ exp) (T) (set A) 


^ P) (T) (set B) 


1 


0.5425 


0.5409 ±0.0044 ±0.0003 


0.5432 ±0.0024 ±0.0002 


2 


0.3551 


0.3619 ±0.0047 ±0.0004 


0.3579 ±0.0036 ±0.0004 


3 


0.1024 


0.0972 ±0.0049 ±0.0004 


0.0991 ±0.0034 ±0.0004 


4 


-0.00055 


-0.0017 ±0.0037 ±0.0004 


-0.0021 ±0.0033 ±0.0003 


5 


-0.4389 


-0.4344 ±0.0050 ±0.0003 


-0.4344 ± 0.0050 ± 0.0003 


6 


0.00067 


0.0037 ±0.0055 ±0.0003 


0.0037 ±0.0055 ±0.0003 


b) Model by Soares [14]: 


i 


bf h) (T) 


Vf xp) {T) (set A) 


b { r P \T) (set B) 


1 


0.3908 


0.3900 ±0.0046 ±0.0003 


0.3955 ±0.0023 ±0.0002 


2 


0.2574 


0.2617 ±0.0049 ±0.0004 


0.2551 ±0.0037 ±0.0004 


3 


0.3518 


0.3483 ±0.0047 ±0.0004 


0.3509 ±0.0037 ±0.0004 


4 


-0.00086 


-0.0083 ±0.0040 ±0.0004 


-0.0017 ±0.0035 ±0.0003 


5 


-0.3171 


-0.3156 ±0.0052 ±0.0003 


-0.3156 ±0.0052 ±0.0003 


6 


0.0011 


0.0008 ±0.0052 ±0.0003 


0.0008 ±0.0052 ±0.0003 






c) Model by Cheng 


[15]: 


i 


bf h) (T) 


bt P) (T) (set A) 


bt P) (T) (set B) 


1 


0.5271 


0.5228 ±0.0045 ±0.0003 


0.5267 ±0.0024 ±0.0001 


2 


0.2928 


0.2980 ±0.0048 ±0.0004 


0.2950 ±0.0036 ±0.0004 


3 


0.1801 


0.1791 ±0.0048 ±0.0004 


0.1778 ±0.0035 ±0.0004 


4 


-0.00066 


-0.0030 ±0.0037 ±0.0003 


-0.0034 ±0.0034 ±0.0003 


5 


-0.3928 


-0.3927 ±0.0051 ±0.0003 


-0.3927 ±0.0051 ±0.0003 


6 


0.00088 


-0.0019 ±0.0054 ±0.0003 


-0.0019 ±0.0054 ±0.0003 
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Table 3: Comparison of the values of observables b\ exp \T), with T' = F s , extracted 
from the Monte Carlo data, with their values bf h ^ (T) corresponding to the model of 
Cheng [15] for initial transversity amplitudes 



% 


bf h) (T) 


bf xp \T) (set A) 


bt P) (T) (set B) 


1 


2.2036 


2.176 ±0.044 ±0.003 


2.206 ±0.026 ±0.001 


2 


1.2242 


1.282 ±0.045 ±0.004 


1.245 ±0.034 ±0.004 


3 


0.9187 


0.917 ±0.045 ±0.004 


0.930 ±0.034 ±0.004 


4 


-0.0073 


-0.099 ±0.036 ±0.003 


-0.094 ±0.032 ±0.003 


5 


-1.6425 


-1.618 ±0.048 ±0.003 


-1.618 ±0.048 ±0.003 


6 


0.0098 


0.067 ±0.050 ±0.003 


0.067 ±0.050 ±0.003 



Table 4: Results of determination of the decay-width parameters (in units (mm/c) -1 ) 
based on extraction of the observables \)f xv ^ from analysis of 100 000 Monte Carlo 
events. The input value of Ar s corresponds to Ar s /T s = —0.15 



Parameter 


Input value 


Measurement (set A) 


Measurement (set B) 


Ar L 


-0.1139 


-0.103 ±0.058 ±0.003 


-0.110 ±0.034 ±0.002 


Ar H 


-0.5696 


-0.478 ±0.137 ±0.012 


-0.554 ±0.101 ±0.012 


r L 


2.4493 


2.444 ±0.029 ±0.002 


2.447 ±0.017 ±0.001 


r H 


2.1076 


2.154 ±0.068 ±0.006 


2.115 ±0.050 ±0.006 


r s 


2.2784 


2.299 ±0.037 ±0.003 


2.281 ±0.027 ±0.003 


Ar s 


-0.3418 


-0.290 ±0.074 ±0.006 


-0.332 ±0.053 ±0.006 
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Table 5: Statistical errors of Ar s extraction (in units [mm/c] -1 ) obtained by applying 
the angular-moments method with set-B (set-A) weighting functions to the Monte 
Carlo data samples with different numbers of events 



Aiyr s 


200 events 


500 events 


10 3 events 


10 4 events 


10 5 events 


-0.03 








0.035 (-) 


0.014 (0.023) 


-0.05 








0.046 (-) 


0.014 (0.022) 


-0.1 






0.11 (-) 


0.046 (0.079) 


0.014 (0.024) 


-0.15 






0.13 (0.19) 


0.045 (0.078) 


0.014 (0.024) 


-0.2 




0.23 (-) 


0.12 (0.18) 


0.048 (0.072) 


0.015(0.026) 


-0.3 


0.21 (-) 


0.23 (-) 


0.18(0.20) 


0.050 (0.083) 


0.016(0.028) 



Table 6: Determination of initial transversity amplitudes and strong-phase difference 
by using the values of observables Vf xp \T) extracted from Monte Carlo data. The 
events sample has been generated for the case of Cheng's model [15] for transversity 
amplitudes and with Ar s = 0.15 Y s . The first errors are statistical while the second 
errors are caused by uncertainties of /S.Y S determination 



Parameter 


Input value 


10 000 events 


100 000 events 


l^o(0)| 2 


0.54 


0.527 ±0.007 ±0.012 


0.5398 ±0.0023 ±0.0011 


|A|,| 2 


0.30 


0.337 ±0.011 ±0.008 


0.3023 ±0.0036 ±0.0006 


|4l| 2 


0.16 


0.136 ±0.010 ±0.020 


0.1579 ±0.0032 ±0.0018 


cos(5 2 - h) 


-1 


-1.021 ±0.044 


-0.9962 ±0.015 
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Figure captions 

Figure 1. Definition of physical angles for description of decays B®(t), B° s (t) — > 
J/^(-> K+K-) in the helicity frame. 

Figure 2. Color suppressed diagrams for decays B® — > J/t/jV ((q, V) G {(s, 0), (d, -ft'*)}). 

Figure 3. Comparison of time and angular distributions for the decays B° s {t) ,~B° s (t) — > 
J/^(-> /i + /u") generated by SIMUB, PYTHIA and QQ packages. For 

time t we use units lmm/c « 3.33 x 1CT 12 sec. 

Figure 4. Dependence of the observables b^T) andS/(T) = &;(T)[l-exp(-rT)]/(TT) 
(i = 1,2,3) on the value of Ar s /r s . The observables have been calculated for the 
case of Cheng's model for transversity amplitudes. 

Figure 5. The flow of data within the package SIMUB. 
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Figure 2 




— i^i^ 

PYTHIA 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 

cos(e K+ ) cos(e^) 




\\ i i i i i i i i i i i I r v •; v r .' j v v i v v v A v v v rv v -y i v , 

2 4 6 0.5 1 1.5 2 2.5 3 

X, rad t, mm/c 



Figure 3 
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Figure 4 



bb_gen bb_frg BB_dec 




Intermediate Ntuple with events 
at parton level in PYJETS format 
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Figure 5 
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